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ABSTRACT 



We report an eclipse mapping study of the intermediate polar DQ Her based 
on time-resolved optical spectroscopy ( AA ~ 3800 — 5000 A) covering 4 eclipses. 
The spectra were sliced into 295 narrow passbands in the continuum and in 
the lines, and the corresponding light curves were analysed to solve for a set of 
monochromatic maps of the disk brightness distribution and for the flux of an 
additional uneclipsed component in each band. Eclipse maps of the He II A4686 
line indicate that an azimuthally- and vertically-extended bright spot at disk 
rim is important source of reprocessing of x-rays from the magnetic poles. The 
disk spectrum is flat with no Balmer or Helium lines in the inner regions, and 
shows double-peaked emission lines in the intermediate and outer disk regions 
while the slope of the continuum becomes progressively redder with increasing 
radius. The inferred disk temperatures are in the range T ~ 13500 — 5000 K 
and can be reasonably well described by a steady-state disk with mass accretion 
rate of M = (2.7 ± 1.0) x 10~ 9 M & yr~ l . A comparison of the radial intensity 
distribution for the Balmer lines reveals a linear correlation between the slope of 
the distribution and the transition energy. The spectrum of the uneclipsed light is 
dominated by Balmer and He I lines in emission (probably from the extended nova 
shell) with narrow absorption cores (likely from a collimated and optically thick 
wind from the accretion disk). The observed narrow and redshifted Call A3934 
absorption line in the total light spectra plus the inverse P-Cygni profiles of the 
Balmer and Hell A4686 emission lines in spectra of the asymmetric component 
indicate radial inflow of gas in the innermost disk regions and are best explained in 
terms of magnetically-controlled accretion inside the white dwarf magnetosphere. 
We infer projected radial inflow velocities of ~ 200 — 500 km s" 1 , significantly 
lower than both the rotational and the free-fall velocities for the corresponding 
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range of radii. A combined net emission He II plus H/3 low- velocity eclipse map 
reveals a twisted dipole emitting pattern near disk center. This is interpreted as 
being the projection of accretion curtains onto the orbital plane at two specific 
spin phases, as a consequence of the selection in velocity provided by the spectral 
eclipse mapping. 

Subject headings: accretion, accretion disks - binaries: eclipsing - novae, cataclysmic 
variables - stars: individual (DQ Her) 
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Introduction 



In many cataclysmic variable stars, the magnetic field of the white dwarf (primary star) 



is sufficiently weak to be neglected. On the other hand, in s ystems 



the magnetic field dominates completely the accretion flow ( IHellier 



ike t 



re AM Her stars, 



20011 ) . Intermediate 



cases can lead to a complex range of phenomena. This occurs in the intermediate polar 
stars, where a magnetic field of moderate strength (B < 10 7 G) allows the combination of 
characteristics of non-magnetic systems (in the outer disk regions) with characteristics of 
highly magnetic systems (in the regions close to the white dwarf). 

In intermediate polars the inner regions of the accretion disk are disrupted by the 
magnetic field of the white dwarf and the infalling plasma is forced to follow the field lines 
down to the star's surface at one or both the magnetic poles, transforming kinetic energy 
into radiation. Part of this energy can be reprocessed in other sites of the binary, leading 
to pulsations at the white dwarf rotation period and/or its beating with the orbital period 



((Warner 



19861 ). 



DQ Her, or Nova Her 1934, is a deeply eclipsing cataclysmic variable (orbital period 
P = 4.6 h) and the prot otype of the subclass of intermediate polars, showing a coherent 71 s 



pulsa tion ( jWalker 



m. 



1956|). Due to the high inclination of the system (i = 86.5°; 



Home et al. 



19931 ) the white dwarf, accretion disk and bright spot are occulted by the secondary star, 
allowing the emission from these different light sources to be se parated and spatially 
resolved by i ndirect imaging techni ques such as eclipse-mapping (IHornd 



tomography ( IMarsh &: Home 



1988 1. 



19851 ) and Doppler 



Here we report the results of an eclipse-mapping study of the spectra and structure of 
the accretion disk of DQ Her based on time-resolved spectroscopy of 4 eclipses collected on 
July 1987. The observations are described in Section [2j Data analysis procedures are given 
in Section |3] In Section H] we investigate the disk spatial structure in the emission lines and 
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in the continuum and we present spatially resolved spectra of the disk, gas stream and of 
the uneclipsed light. The results are discussed in Section [5] and summarized in Section O 
In the Appendix |A] we discuss the ability to reconstruct the full-width-half-maximum 
distribution of emission lines with spectral mapping techniques. 

2. Observations 



The data consist of 3272 spectra obtained with the Double Spectrograph at the Hale 
5 m telescope, using the 2D-FRUTTI detector on the nights of 1987 July 3 — 4. The 



runs cover the ec 



Africano fc Olson 



i pse cy cles 62105, 62106, 62110 and 62111, according to ephemeris of 



( 1l98ll ). The integration time for each spectrum was 10 s, and the usa ble 



spectral coverage is AA ~ 3800 — 5000 A. The reader is referred to 



Martell et al. 



( 19951 ) for 



a detailed description of the observations and of the data reduction procedures. 

The data were correcte d to the rest fram e of the white dwarf by removing its orbital 



velocity K\ = 140 kms 1 ( JHorne et al.lll993l ) from the spectra, and then binned to a 



resolution of 130 kms -1 per bin (2 A pixel -1 at 4500 A). In the present work we analyse 
only the data around eclipse, in the phase range —0.1 to +0.1, comprising a subset of 1073 
spectra from the original data set. 

Fig. [Q shows average spectra prior to (phase range —0.09 to —0.08 cycle), after (+0.10 
to +0.11) and at mid-eclipse (—0.02 to +0.02). The spectra show Balmer emission lines, a 
strong He II 4686 line, weak He I and He II lines and a few metal absorption lines on top of 
a flat continuum. The emission lines show a clear double-peaked profile at mid-eclipse. The 
Balmer lines are more prominent at mid-eclipse than outside eclipse, indicating that their 
emission is less concentrated towards disk center than that of the continuum. 
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Fig. 1. — Average ingress (dark grey, phase range —0.09 to —0.08 cycle), egress (light grey, 
+0.10 to +0.11 cycle) and mid-eclipse (black, —0.02 to +0.02 cycle) spectra of DQ Her, from 
the combined eclipse cycles 62105, 62106, 62110 and 62111 (see text). Vertical dotted lines 
indicate the emission/absorption lines. Horizontal ticks mark the 295 passbands used to 
extract light curves. For a better visualization the ingress and egress spectra are vertically 
displaced by —1.4 and 1.7m Jy, respectively. 
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3. Data analysis 



3.1. Light-curve construction 



The spectra were divided into 295 narrow passbands; 72 continuum passbands, typically 
8 A wide, and 223 velocity-resolved bands of 200 km s~ l for the 22 lines marked in Fig. [TJ 
In defining the velocity-re solved line bins, w e neglected the small systemic velocity of the 



object (7 = 5 ± 5 km s 1 . 



Home et al. 



19931 ). For the passbands including emission lines, 



the light curve comprises the total flux at the corresponding bin, with no subtraction of a 
possible continuum contribution. 

Light curves were extracted for each band, combining the average flux of the four data 
sets on the corresponding wavelength range in phase bins of 0.0025 cycles. The error bars 
are taken as the standard deviation with respect to the average flux at each phase bin. T he 



result ing light curves were phase-folded according to the ephemeris of 



Africano fc Olson 



(1981). 



T mid = HJD 2434954.943899(2) + 0.193620897(5) E. 



;i) 



We remark that, because the 71 s oscillation was not removed from the light curves, the 
resulting eclipse maps contain the combined light distribution of the steady, unpulsed light 
and a time averaged distribution of th e pulsed light. Ney ertheless, since the 71 s accounts 



for, at most, 2 — 3 % of the total flux (IMartell et al- 



ii}^ , the influence of the pulsed light 



on the overall disk brightness distribution is small. Light curves and eclipse maps for the 
H/3 and He II 4686 lines are shown in Figs. [2] and [3l respectively. Contour curves in these 
figures enclose the regions of the eclipse maps at and above the 15 a level of statistical 
significance. 



3.2. Eclipse-mapping 



We used maximum-entropy eclipse- mapping techniques ( JHornelll985l ; iBaptista &: Steiner 



19931 ) to solve for a map of the disk brightness distribution and for th e flux of an addi tional 



Baptistal (j200lh for a 



uneclipsed component in each spectral band. The reader is referred to 
review on the eclipse-mapping method. 

As our eclipse map we adopted a flat grid of 65 x 65 pixels centered on the primary 
star with side 2 Rl%, where Rn is the distance from the disk center to the inner Lagrangia n 



point (Li). In this work we adopted a Roche lobe size of Rli = 0.766 R & (IPetterson 



1980|) 



The eclipse geometry is defined by the mass ratio q, which sets the shape and the relative 
size of the Roche lobes, and the inclination i, which determines the shape and extension of 
the shadow of the secondary star projected onto the orbital plane. We adopted the values 



by 



Home et al. 



( 119931 ). q = 0.66 ±0.04 and i = 86.5° ±1.6°. This combination of parameters 



ensures that the white dwarf is at the center of the map. 

Before applying eclipse mapping techniques to a high inclination system such as 
DQ Her, it is important to address the question whether the disk rim can hide parts 
of the accretion disk from view (thereby leading to unreliable results). In particular, a 
vertically-extended, thick disk rim could lead to an important self-obscuration effect on the 
disk side closest to the observer if the disk half-opening angle is (3 > (90° — i). An estimat e 



of the disk half-opening angle in DQ Her comes from the recent work of 



Saito fc Baptista 



( 120091 ). They applied eclipse mapping techniques to trace the changes in surface brightness 
distribution along the 71 s pulsation cycle and found that a vertically-extended bright spot 
at disk rim is the main site of emission of the 71 s optical oscillations. The upper side of the 
bright spot is illuminated by one of the magnetic poles once per cycle, while the lower side of 
the bright spot is illuminated by the other pole roughly half a cycle later. From the difference 
in the projected position of the upper and lower bright spot sides in their eclipse maps, 
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4 



they inferred a vertical thickness of H = (7.4 ± 0.4) x 10 3 R & [tan(i)/16.35] : HI, where the 
term within brackets is unity for i = 86.5°. Combined with a bright spot radial distance of 
R BS = (0.44 ±0.02) R & (see SectionE]), this leads to = arctan{H/2R BS ) = 0.5° - 0.002° 
for i = 86.5° — 89°. The lowest range of values for (3 is hard to be interpreted by the 
steady-state disks model, where flared disks are expected, with the vertical tickness H 
proportional to the mass accretion rate M and increasing with radius with by factor 
H oc r 9//8 . The inferred value of (3 decreases with increasing inclination, and the inequality 
f3 < (90° — i) holds for any reasonable inclination (i > 70°). Hence, self-obscuration by 
the disk rim is not an issue in DQ Her. The fact that the bright illumination pattern that 
rotates with the 71 s pulsation period in t he inner disk regions re mains visible when it is 



passing in the disk side near the observer (ISaito fc Baptista 



20091 ) yields additional support 



for the above statement. We are therefore confident that no distortion effects caused by 
self-obscuration by the disk rim arise when applying eclipse mapping techniques to the 
DQ Her light curves. 

For the r econs tructions we used the default of limited azimuthal smearing of 



Rutten et al. 



(119921 ). which is better suited for reco vering asymmetric st ructures than 



the original default of full azimuthal smearing (c.f., 



Baptista et al. 



19961 ) . Simulations 



showing the ability of the eclipse mapping method to reconstruct asymmetric brightness 



distributions (even under unfav orable condition s of l ow S/N and incomp 



are presented and discussed by 



Baptistal (I200l[ ) and 



e te pha se coverage) 



Baptista fc Catalan! (J2Q0lh . We used a 



radial blur width Ar = 0.0185-Rli and an azimuthal blur width A# = 20°. 

The statist ical uncertainties in the eclipse maps were estimated with a Monte Carlo 
procedure (e.g., iRutten et al.lll992l ). For each narrow band light curve, a set of 20 artificial 



Tn 



Saito fc Baptistal (120091 ) this value is mistakenly quoted as the bright spot vertical 



half-thickness. 
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light curves was generated in which the data points were independently and randomly varied 
according to a Gaussian distribution with standard deviation equal to the uncertainty at 
that point. The light curves were fitted with the eclipse-mapping code to produce a set of 
randomized eclipse maps. These were combined to produce an average map and a map of 
the residuals with respect to the average, which yields the statistical uncertainty at each 
pixel. The uncertainties obtained with this procedure were used to estimate the errors in 
the derived radial intensity and temperature distributions as well as in the spatially resolved 
spectra. We also divided each eclipse map by the map of the residuals with respect to the 



average to produce maps of the inverse of the relative errors, or signal-to-noise ratio maps 
( jBaptista fc Bortolettol 



2004]). The S/N maps are overplotted on the corresponding eclipse 
maps in Figs. [2] and |3] as a contour line for S/N = 15. All structures in the eclipse maps are 
significant at a confidence level equal to or above 15a. 

4. Results 

4.1. Disk structures 

Fig. [2] shows eclipse maps of H/3 for a set of 5 velocity bins. Similarly, Fig. [3] shows 
velocity-resolved eclipse maps for the Hell 4686 line. The velocity-resolved H/3 (and, to 
a less extent, the Hell 4686 line) light curves show the expected behavior for the eclipse 
of gas rotating in the prograde sense (rotational disturbance), with the blue side being 
eclipsed earlier than the red side. The eclipse maps in symmetric velocity bins do not show 
the reflection symmetry around the line joining both stars expected for the emission from a 
symmetric disk around the white dwarf. 

The clear egress shoulder in the light curves reveals an azimuthally-extended asymmetry 
in the outer parts of the disk extending ~ 90° ahead of the intersection of the stream 
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Fig. 2. — Data and model light curves and eclipse maps for H/3 (A4861) at five velocity bins 
(v = ±800, ±400, kms' 1 , Av = 200 kms~ 1 ) . The left panels show the data (dots with 
error bars) and model (solid lines) light curves. The middle panels show the corresponding 
eclipse maps for the total emission (line ± continuum) in a logarithmic greyscale. Brighter 
regions are indicated in black; fainter regions in white. A cross marks the center of the 
disk; dotted lines show the Roche lobe, the gas stream trajectory and a circle of radius 
R-ss = 0.57 Rl\ corresponding to the bright spot position; the secondary is to the right of 
each map and the stars rotate counter-clockwise. Contour curves enclose the regions at and 
above the 15 a level of statistical significance. The right panels show eclipse maps of the 
net line emission in each case; the notation is similar to that of the middle panels. The 
horizontal bars indicate the logarithmic intensity level of the greyscale in each column. 
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Fig. 3. — Velocity-resolved light curves (left) and eclipse maps (middle and right) for Hell 
4686. The notation is similar to that of Fig. [2j 
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trajectory with the disk rim. This azimuthally-extended bright spot indicates that the shock 
of the inflowing gas with the outer disk is radiatively inefficient, with the emission pattern 
spreading along the direction of gas rotation. Moreover, if the bright spot cannot cool 
efficiently it should also expand in the vertical direction, leading to a vertically-extended 
bulge in the outer disk that becomes a natural site for re processing of x-ray radia tion from 



the rotating magnetic poles - in line with the findings of ISaito fc Baptistal (120091 ). From 



the radial position of the maximum intensity of the bright spot, we estimate an outer disk 
radius of R BS = 0.57 ± 0.03 R Ll . 

In order to emphasize the structures observed in the line maps, we calculated net line 
emission maps. These were obtained by combining continuum eclipse maps on the short 
and long wavelength sides of the target emission line and by subtracting the derived average 
continuum map from each of the velocity resolved line maps. Because of our choice to 
display the maps in a logarithmic grayscale (that only handles positive intensities), regions 
where the resulting intensities are negative (indicating line absorption) were suppressed 
from the net line emission maps shown in Figs. |2] and [3j The Balmer lines net emission 
maps show a dearth of emission at the bright spot position, indicating that these lines are 
in absorption at the infalling stream impact site (see Section H~2j) . In contrast, the Hell 
4686 line appears clearly in emission at the bright spot position (Fig. EJ right-hand panels). 
The emission is stronger in the positive velocity maps but weak blueshifted emission is also 
seen in the —400 and —800 km s _1 velocity maps. The latter is likely the result of gas 
splashing away from the accretion disk (with a velocity component towards the observer) 
after the collision of the infalling gas stream with the disk rim. The redshifted emission 
reflects the fact that disk gas at the bright spot position is mostly moving towards the far 
side of the disk and, therefore, away from the observer. 

Because the temperatures in the outer disk regions (T ~ 6000 — 7000 K, see Section |4~5|) 
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are not enough to power the high-excitation He II 4686 line, the strong He II emission from 
this region is an indication that the bright spot must be an important source of reprocessing 
of x-ray radiation from the magnetic poles/curtains - again in line with the findings of 



results of 



Saito fc Baptist a (2009J). This is consistent with and leads to a natural explanation for the 



Patterson et al. 



(Il978[ ) -- who found that the amplitude of the 71 s pulsations is 



variab le along the orbit and reaches maximum around phase <fi ~ 0.25 — and of lO'Donoghue 



( 119851 ) - - who correctly advanced the conclusion that the orbital behavior of the 71 s 
pulsations amplitude requires that the x-ray reprocessing site be a disk rim with variable 
thickness and maximum visibility at phases 0.2 — 0.3. This coincides with the phase range 
where the upper side of the azimuthally- and vertically-extended bright spot is best seen. 

The zero- velocity H/3 net line emission map shows a peculiar "heart" -shaped 
distribution close to disk center (R < 0.3 Rli) with major axis aligned at an angle of ~ 150° 
with respect to the line joining both stars (Fig. [2]). Morphologically similar structures 
are also seen in the ±200 and ±400 kms~ l velocity maps, with the positive velocity 
maps showing the receding part of the structure (i.e., the upper hemisphere of the eclipse 
maps in Figs. [2] and |3]) and the negative velocity maps showing the leading part of the 
structure (i.e., the lower hemisphere of the eclipse maps in Figs. [2] and [3]). An asymmetric 
structure resembling a dipole pattern and aligned at a similar orientation can be seen in 
the zero- velocity Hell 4686 net emission map (Fig. [3]). We will return to this point in 
Section 15.21 



4.2. Spatially resolved spectra 

Motivated by the distinct emission observed at the bright spot and gas stream regions 
and in order to avoid contamination by bright spot emission that could affect the spectrum 
of the near side of the disk (the hemisphere closest to the secondary star in the eclipse 
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maps of Figs. [2] and [3]), we divided the eclipse maps in two distinct regions: "disk" and 
"gas stream". We define the "disk" as the disk section between azimuths 90° and 270°, 
and "gas stream" as the region between azimuths 0° and 90° (see Fig. H]). We further 
defined annular sections of width 0.1 Rli to extract spatially-resolved spectra as a function 
of radius. In order to isolate the contribution of possible asymmetric brightness sources 
in the eclipse maps (e.g., the bright spot, gas stream or accretion curtains) we separated 
the total emission into its symmetric and asymmetric disk emission components. The 
symmetric component is obtained by slicing the disk into a set of radial bins, computing 
the average of the lower half of the intensities in each bin, and fitting a smooth spline 
function to the resulting set. The spline fitted intensity in each annular section is taken 
as the symmetric component. The asymmetric component is obtained by subtracting the 
symmetric component from the original eclipse map. This procedure essentially preserves 
the baseline of the radial profile, removing all azimuthal structure from the symmetric 
component. The statistical uncertainties affecting the fitted intensities are estimated with a 
Monte Carlo procedure (Section 13.21) . 

Fig. [5] shows spatially resolved spectra for the symmetric disk component and for the 
"gas stream" region (total emission). In the inner regions (R < 0.2 Rli) the disk spectrum 
shows a flat and featureless continuum with no Balmer or Hell lines. The disk spectrum 
becomes progressively redder and fainter with increasing radius, indicating the existence of 
a radial temperature gradient. The Balmer and Hell lines start to appear in emission at 
intermediate disk regions (R ~ 0.3 Rli) and become stronger with increasing radii. The 
spectrum of the "gas stream" is systematically bluer and brighter than the disk spectrum at 
the same radius and the difference increases with radius. In contrast to the observer in the 
disk spectra, in the outer regions (R > 0.5 -Rli) the lines appear in absorption, indicating 
that the emitting gas is optically thick. The P Cygni profile of the Balmer lines in the "gas 
stream" spectrum framing the bright spot, at R = (0.5 — 0.6) Rli, is an indication that 
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part of the colliding gas stream deflects back towards the observer, in agreement with the 
inference previously drawn from the He II A4686 net emission line maps. 

The lack of Balmer lines in the innermost disk regions may be interpreted in a scenario 



of magnetically-controlled accretion for radii sma 



for DQ Her a white dwarf mass of M\ = 0.6 M & ( iZhang et al. 



ler than the Alfve n radius. Assuming 



19951 ). an accretion rate of 



M = 2.7 x 10 9 M & yr 1 (Section 14. 5p and a white dw arf magnetic mo ment in the range 



Patterson 



1994 ). we estimate an 



fi = 10 32 — 10 33 G cm 3 (typical for intermediate polars, 
Alfven radius for DQ Her in the range Ra — (0.05 — 0.25) Rli- It is expected that inside 
this radius the disk will be disrupted, with the gas being channelled by the magnetic field 
lines onto accretion curtains down to the magnetic poles at the white dwarf surface. This 
produces magneto-bremsstrahlung (free-free) radiation by the optically thin gas, which leads 
to a spectrum basically flat for the spectral range between the self-absorption limit and a 
thermal-cutoff. Hence, the spectrum observed in the inner most disk regions (i.e., flat an d 



featu reless) could plausibly be interpreted in this context (jKraft 



1988|). 



1959; 



Wickramasinghe 



A relatively narrow (FWHM < 400 kms^ 1 ) and redshifted (by ~ 200 kms^ 1 ) Call 
A3934 absorption line can be seen both in the disk and in the gas stream spectra up to 
R ~ 0.3 Rli- (There are hints that a similarly narrow and redshifted Call A3968 absorption 
line is also present in the disk and gas stream spectra at the same radial range, although 
we assign a low weight to this statement as this line is heavily blended with the broad 
He emission line.) This corresponds to the Call s-wave seen in absorption in the trailed 



Martell et al. 



(119951 ) associated to the accretion 



spectrogram of the same data set and which 
disk. 

Could the Call line be the result of absorption of disk radiation by a cool and 
vertically-extended disk rim along the line of sight? The evidences say no. If this was the 
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case, the line strength would increase with radius for spectra of the near side of the disk 
since the obscuration effect increases as one moves towards the disk rim in the near side. 
On the contrary, the observed Call line strength decreases with increasing radius in the 
spectra of the "gas stream" (Fig. [5]) as well as in the symmetric component of spectra 
extracted for the near side of the disk (not shown here). 

The temperatures inferred for the inner disk regions (T > 10 4 K, see Section I4.5[) 
are too high to allow for significant Call emission. Thus, the Call absorption should 
be produced in intervening cooler gas along the line of sight. Interesting inferences can 
be drawn from the dynamics of this intervening gas. Gas rotating in Keplerian orbits at 
this distance from the white dwarf should lead to broad (FWHM ~ 1000 — 3000 kms^ 1 ) 
and double-peaked absorption lines with no blueshift/redshift. The observed narrow, 
single-peaked and redshifted Call line indicates that the gas dynamics in the absorbing 
region is dominated by radial motion (in opposition to azimuthal, Keplerian motion). 
Depending on where the absorption occurs along the line of sight, we may be witnessing 
gas outflow (if absorption occurs closer to the orbital plane, at the far side of the disk) or 
gas inflow (if absorption occurs closer to the observer along the line of sight, at the near 
side of the disk) . 



In the radial outflow scenario DQ Her would be acting as a propeller (e.g., 



Wynn et al. 



19971 ). expelling gas from inside the magnetosphere towards larger radii. The corotation 
radius for DQ Her (R co = 0.031 R & = 0.041 Rli) is comparable to (and even smaller than) 
our lower estimate of the Alfven radius, suggesting that plas ma blobs that attach to field 

19951 Chapter 7). 



lines at R co < R < Ra should be centrifugally accelerated ( jWarner 
Accordingly, one would expect to see blueshifted Ca II absorption in spectra of the near 
side of the disk in account of gas expelled towards the observer. However, contrary to 
this expectation, the Call absorption line is also redshifted both in the "gas stream" 
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spectra and in the symmetric component of spectra extracted for the near side of the disk. 
Additional problems with this scenario include the difficulties to explain why the outflow 
velocity would be so low (v out — 200 km s" 1 ), why it does not change with radius, and why 
no Call absorption is seen in the outer disk regions. 

The radial inflow alternative leads to a magnetically-controlled accretion scenario, in 
which the Ca II absorption line is produced when relatively cool regions of the accretion 
curtain are seen projected against hot gas in the background. Inside the magnetosphere 
(R < Ra) two effects contribute to broaden the emission/absorption lines: (i) the gas 
corotates with the white dwarf magnetic field, the spin velocity of which increases with 
radius, v $ = R Co> sp in ; (ii) at the same time, gas moves inwards along field lines with a 
radial velocity v r < Vff oc Rr x l 2 . Because DQ Her is a fast rotator, large azimuthal 
velocities (y^ ~ 2300 — 12000 kms~ l ) are expected for the range of radii of interest, 
R = 0.05 — 0.25 Rli, arid lines produced in this region should be quite broad. Then, why is 
the Call line so narrow? 

At this point it is worth noting that because an eclipse map yields a snapshot of the 
eclipsed brightness distribution over a time scale much longer than the 71 s pulsation period, 
the contribution of the rapidly rotating magnetosphere to our spatially-resolved spectra is 
a time averaged spectru m over all orientations ( or, pulse phases). Adopting the oblique 



dipole accretor model of (jSaito &: Baptista 



20091 . see their Fig. 3), significant absorption 



should occur around spin phase spin = 0.5 — when the upper accretion curtain covers the 
largest portion of the hot inner magnetosphere and/or the (also hot) innermost regions of 
the far side of the accretion disk. However, at this orientation the spin component of the 
accretion curtain velocity field is mostly in the plane of the sky and the projected velocity 
along the line of sight is essencially v r . Note that, at quadrature (0 sp in = 0.25 and 0.75), 
the absorption of radiation by the upper accretion curtain is not only significantly reduced 
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but is also smeared in wavelength because of the contribution of a large spin velocity 
component to the line profile. Thus, the observed absorption line profile is dominated by the 
spectral contribution of the accretion curtains at sp i n = 0.5. The magnetically-controled 
accretion scenario implies radial inflow velocities of v r ~ 200 km s _1 , significantly lower 
than both the rotational and the free-fall velocit i es for that range of radii. This is line 



with inferences drawn by IWu fc Wickramasinghd (Il99ll ). who pointed out that, for rapidly 



rotating intermediate polars, the accretion shock may occur almost at Ra and the infall 
from Ra to the white dwarf surface is unlikely to have a significant radial component. 

4.3. The emission lines 

In this section we analyse the radial behavior of Balmer (H/3, H7 and US) as well as 
the He II 4686 lines. Fig. [6] shows the radial intensity distribution for the lines and for the 
adjacent continuum (top), the radial intensity distribution for the net line emission (second 
row from top), the radial run of the equivalent width (EW, third row from top) and the 
radial run of the full width at half- maximum (FWHM, bottom). The diagrams were derived 
from the symmetric disk component (Section 14.21) . The performance of the eclipse mapping 
technique to reconstruct FWHM distributions is discussed in the Appendix [A] 

The intensities in the line maps are higher than in the adjacent continuum at 
intermediate and outer disk regions, indicating that the lines appear in emission there. The 
Balmer lines disappear in the continuum in the innermost disk regions (R < 0.2 Rli), and 
the corresponding net line emission becomes unreliable or meaningless. For H7 and B.S the 
line intensity is lower than the continuum level at disk center indicating that these lines are 
in absorption there. These regions are not plotted in the corresponding net emission panels. 

The net line emission decreases in strength with increasing radius in the intermediate 
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and outer disk regions. The radi al dependence has a particular slope for each line, indicating 



that the empirical law I oc r 5 (jMarsh et al 



19901 ) is not adequate to describe the behavior 



of the emission lines in DQ Her. Moreover, we found that, for the Balmer lines, the radial 
dependence is linearly correlated with the energy required to excite the line. This result 
will be discussed in Section 15.11 

The EW of the Balmer lines increase with radius with maxima of EW ~ 100 — 400 A at 
R = 0.5 Rli- In contrast, the Hell 4686 line shows an EW ~ 20 A, roughly independent of 
radius. In all cases, the EW becomes negligible (or negative) at the innermost disk regions 
because the continuum intensities reach (or exceed) the line intensities. 

In the intermediate and outer disk regions (R > 0.3 -Rli) the slope of the observed 
FWHM distribution of the Balmer lines is fairly consistent with the expected vk oc R~ 1 ^ 2 law 
for gas rotating in Keplerian orbits. However, the observed FWHM values are systematically 
lower than the predicted for a disk around a white dwarf with M\ = 0.6 M Q (the dashed 
line in the lowermost panels of Fig. [6]) by ~ 30-40%, suggesting that either the adopted 
white dwarf mass is overestimated or that the emitting gas has sub-Keplerian velocities. No 
reasonable error in the estimated inclination could account for this discrepancy. Moreover, 
in order to reconcile the Keplerian velocity distribution with the observed FWHM values, 
an unrealistically low-mass white dwarf (< 0.35 M & ) would be needed. It is therefore likely 
that the disk gas has sub-Keplerian velocities. The He II 4686 FWHM is hard to measure 
and becomes unreliable for R > 0.3 Rli, because the net line intensity is decreasing sharply 
with radius and the EW is relatively small. 

Fig. [7] shows the velocity profile for the total and the asymmetric component of the 
emission lines. In the total emission spectra the lines appear in emission in the intermediate 
and outer regions (R > 0.2 Rn) becoming progressively sharper (stronger and narrower) 
with increasing radius. The Balmer and the Hell 4686 lines show a clear double-peaked 
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profile with the blue peak generally stronger than the red peak in the Balmer lines. The 
Hell line changes from a double-peaked profile at small radii to a single-peaked profile in 
the outer disk regions. 

The differences in behavior between the Balmer lines and the Hell 4686 line has 



already b een observed in 



is absent ( IChanan et al 



previous analyses o 



19781 ) or negligible ( IMartell et al. 



DQ Her. The 71 s pulsation, evident in He II, 



19951 ) in the Balmer lines. This 



result reflects the different ionization/excitation conditions between hydrogen and heavier 
ions. 

The spectra of the asymmetric component evince inverse P-Cygni profiles for all lines, 
with a 'red' absorption component on top of a 'blue' emission peak. The strength of the 
absorption component decreases with increasing radius and disappears for R > OAR^i. 
The Hf3 absorption component is narrow and slightly blueshifted (by ~ — 200 km s^ 1 ) in 
the innermost regions, and progressively broadens and moves towards positive velocities 
with increasing radius, reaching trough velocities of ~ +500 km s -1 at R = (0.3 — 0.4) Rli- 
The absorption component is always redshifted for the other lines, also with velocities of 
~ +500 ferns -1 . The emission components are centered at ~ —500 /cms -1 . The blueshifted 
emission seen in spectra of the asymmetric component of the far side of the disk indicates 
that emitting gas of a non axi-symmetric brightness source is falling radially towards the 
white dwarf at disk center. Consistently, absorbing gas along the line of sight (possibly at 
the near side of the disk) is also falling radially towards the white dwarf at disk center with 
similar radial inflow velocities (v r ~ 500 ferns -1 ). These evidences are in agreement with 
those inferred from the behavior of the Ca II A3934 absorption line; they give additional 
support to the magnetically-controlled accretion scenario put forward in the preceeding 
section. 
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4.4. The uneclipsed component 

At the high inclination of DQ Her, there is no region in the far side of the disk which 
is never eclipsed and, therefore, it is not possible to dump extra flux at mid-eclipse without 
affecting the shape of the model eclipse light curve. Any extra flux at mid-eclipse is 
unambiguously attributed to 'uneclipsed light', i.e., the fraction of the total light which 
cannot be included in the eclipse map. Thus, the uneclipsed component in DQ Her arises 
from light outside of its orbital plane (or from the outward-facing hemisphere of the 
eclipsing secondary star). 

Fig. [8] shows the uneclipsed component and its fractional contribution with respect 
to the average out-of-eclipse level. The uneclipsed spectrum is dominated by the Balmer 
lines, which appear in emission with an apparent double-peaked profile. H/3 is the most 
prominent line, reaching ~ 25 % of the total emission at its wavelength. The He I A4921 
line is also strong in the spectrum, at a level similar to that of H<5 but with a wider line 
profile. In contrast to the observed in the average spectra (in- and out-of-eclipse), where 
the continuum is almost flat, the continuum in the uneclipsed component rises towards 
longer wavelengths, reaching the same intensity level of some lines (e.g., the He line). 
A close inspection indicate that the uneclipsed Balmer lines have a shallow and narrow 
central absorption core (FWHM ~ 200 km s^ 1 ) on top of a broader emission component 
(FWHM ~ 800 /cms" 1 ). The absorption core of the He I A4921 line is deeper and wider, 
with a FWHM ~ 600 km s -1 . 

It is unlikely that the uneclipsed emission lines are from the secondary star. Although 
trailed spectra of Balmer lines show an emission component with the secondary star motion, 
the emission comes from the inward-facing hemisphere of the secondary star — indicating 
it is caused by irradiation effects -- and is apparently restricted to orbital phases around 



0.5 — suggesting that the emitting region is not visible during eclipse ( JMartell et al 
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19951 ). As the inward- facing hemisphere of the secondary star is not seen around eclipse, 
it should not contribute to the uneclipsed light. Furthermore, the width of the uneclipsed 
lines is much larger than expected for a r otationally broade ned line from the secondary star 



in DQ Her (Vrotsinz = 115 ± 18 km s l , 



Home et al. 



19931). 



The uneclipsed emission lines are probably from the nova shell. We note that nebular 
emission from the extended nova shell was included in the slit together with the light 



from the binary, and its contribution was not removed from the spectra ( IMartell et al. 



19951 ). Accordingly, the FWHM of the uneclipsed emission lin es is consistent wit h the 



Vaytet et al 



20071 ). The 



expansion velocity of the nova shell (v = 370 ± 14 /cms -1 , see 

narrow absorption component possibly arises in a collimated and optically thick wind from 

the accretion disk. Evidences for the presence of a collimate d disk wind in DQ He r come 



from the eclipse behavior of the CIV A1550 resonan ce line ( lEracleous et al 



study of scattered x-ray emission (IMukai et al. 



19981 ). the 



20031 ). and from t he ellipsoidal 



shape of the 



2007m . The 



nova shell and the ablation of gas clumps around the shell poles (jVaytet et al. 
secondary star may be responsible for the rising continuum at the red end of the spectrum, 
although the short wavelength range where the continuum appears prevents a meaningful 
attempt to infer its spectral type. 



4.5. The temperatures of the disk 

Here we fit isothermal, pure hidrogen atmosphere models (HI), to the spatially resolved 
spectra at each disk annular region (Fig. [5]) in order to derive the radial dependence of 
temperature and column density in the disk (Fig. |9]). We used the symmetric component of 
the emission and we masked the regions of the lines, since the models only account for the 
continuum emission produced by bound-free and free-free H I transitions. The fit depends 
on the solid angle 9 2 assumed for the pixel in the eclipse maps, which by its turn depends 
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on the distance and the inclination with which the pixel is seen by an observer at the Earth. 
The effective inclination of a pixel in the far side of the accretion disk is W/ = i — 0, where 



is the dis k half-opening ang 



and 525 pc ( iVaytet et al. 



e. We obtained fits for distances of 400 pc fjHorne et al. 



19931 ) 



20071 ). and half-opening angles of = 0°, 2°, and 3.5°. 



Fig. shows the resulting radial dependence of the temperature and column density 
for the case d = 400 pc and = 2° (solid line with filled triangles) and d = 525 pc and 
= 2° (open circles). Also shown are the results for two extreme cases of low (larger 
distance and lower 0) and high (lower distance and larger 0) 9 2 . The \ 2 °f the fit increases 
perceptibly for d = 525 pc, because the high temperatures required to match the average 
flux level of the spectra can no longer reproduce the red slope of the observed spectra. 
Although the different choices of 9 2 lead to systematic differences in temperature and 
column density, these differences are of the order of the uncertainties affecting the derived 
quantities and do not alter the inferences that can be drawn from Fig. [9] Our results 
(for d = 400 pc and a = 2°) indicate that the surface density increases with increasing 
radius, showing values of E ~ 4 x 10~ 5 gcm~ 2 at R = 0.15 Rli and S~3x 10 2 gcm~ 2 at 
R = 0.65 Rli- The steady-state disk models predict that the surface density is proportional 
to a~ 4 / 5 -R~ 3//4 , where a is the viscosity parameter. Thus, our values for £ are consistent 
with steady-state model only if the viscosity in the disk decreases with increasing radius 
in order to compensate the factor oc i?~ 3 ' 4 . The brightness temperature ranges between 
T ~ 13500 K at R = 0.15 R L1 and T ~ 5000 isT at R = 0.65 R L1 , with a radial dependence 
reasonably well described by a steady-state optically thick disk model with mass accretion 



rate of M = (2.7 ± 1.0) x 10" 9 M yr~ l The inferred 



agreement with th e values derived by (jChanan et al. 



( iZhang et al. 



mass accretion rate is in reasonable 
M = 10~ 9 M yr" 1 ) and 



1978 



199.4 M = 3.4 x 10~ 9 M Q yr- 1 ) 



Temperatures in excess of 8000 K are required in order to produce the strong lines 
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with shallow Balmer decrement observed in the intermediate and outer disk regions. This 
is in contrast with the low temperatures (T < 5500 K) inferred from the slope of the 
red continuum emission, suggesting the existence of a vertical temperature gradient in 
the DQ Her disk gas, with a hot, optically thin chromospheric layer (responsible for the 
emission lines) on top a cool, opaque and dense disk photosphere (responsible for the red 
continuum) . The chromosphere is likely heated by irradiation of x-ray photons from the 
magnetic poles and accretion curtains, leading to high-excitation lines. This scenario is 
beyond LTE simple models. 



Smakl (Il98ll ) remarks that in DQ Her the continuum emission seems to be a function 
of the accretion disk regime while the line emission appears to be correlated with the 
luminosity or surface temperature of the primary star. In this context the lines are formed 
by reprocessing of the energy radiated from the central source, probably the accreting 
material falling over the magnetic poles of the white dwarf. T his central emission in DQ 



Her can be detected in soft x-ray emission (IMukai et al 



20031 ) 



We note that the temperatures in the outer disk regions ( R > 0.3 fin) a re systematically 



below the critical temperature for disk instabilities to set in ( iWarner 



19951 ) . suggesting that 



the DQ Her disk could show dwarf nova-type ou tbursts. This holds even if one assumes the 



Zhang et al. 



(119951 ) (illustrated by the steady-state disk model 



higher mass accretion rate of . 

with M = 10~ 8 ' 5 in Fig. [9]). There is no record of dwarf nova outbursts in this binary. 

We end this section raising a puzzlingly different scenario from the radial temperature 
distribution of Fig. O Given the evidences presented here of magnetically-controlled 
accretion for R < 0.3 Rli, it is reasonable to assume that a bona fide accretion disk exists 
only for distances outside this range, and that the brake in slope of the temperature and 
column density distributions at R ~ 0.3 Rli reflects the change from the azimuthal flow 
of gas in the outer accretion disk to a magnetically-controlled radial flow of gas over the 
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accretion curtains in the inner regions. In this case, the outer accretion disk s hows a flat 
temperatur e distribution resembling those found in quiescent dwarf novae (e.g., 



Wood et al. 



1986 



19921 ). The inferred low temperatures (T ~ 6000 f^) and high surface densities 
(S ~ 10 4 g cm~ 2 ) suggest that the disk viscosity is low. This is in marked contrast to the 
common view that the old nova DQ Her is a nova-like variable with a hot, viscous accretion 
disk. 



5. Discussion 



5.1. An energy dependency of the emission-line radial intensity distribution 

slope 

In Section fl~3l we found that the slope of the radial line intensity distributions correlates 
with the transition energy of the Balmer lines. In order to extend our analysis we extracted 
(using the same procedures described in Section 14. 3p the radial intensity distribut ion for 



He, the next line in the Ba lmer series. We then fitted a power-law I oc r a (e.g., 



1981 



Home fc Marsh 



Smak 



19861 ) for each radial distribution in order to find the best-fit a index 
for each line. The innermost disk regions (R < 0.2 Rli) are excluded from the fit since 
the continuum reaches the emission level of the lines in the regions inside this radius. The 
results are listed in Table 1 and plotted in Fig. [TUJ 

The a index increases with the transition energy for the Balmer lines, i.e., the emission 
at lines with larger transition energy (e.g., He) are more concentrated at the disk center 
than lines of lower energy (e.g., H/3). Moreover, our results reveal that there is a strong 
linear correlation between a and the transition energy fo r the Balmer lines (see Fig. [TOl) . 



This result is in contrast with the empirical I oc r law (iMarsh et al. 



199(3 ) usually evoked 



in works about radial distribution of emission lines in accretion disks. The He II 4686 line 
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Fig. 4. — Schematic diagram showing the regions defined as "disk" and "gas stream" . Dashed 
lines mark the annular regions of width 0.1 Rli used to extract spatially resolved spectra. 
Dotted lines show the projection of the primary Roche lobe onto the orbital plane and the 
gas stream trajectory. Azimuths are measured with respect to the line joining both stars 
and increase counter-clockwise. Four reference azimuths are labeled in the figure. 



Table 1: Index a for the Balmer series and Hell lines. 



Line 



H/3 
H 7 
R5 
He 



E[eV] a {I oc r~ a ) 



12.76 1.20 ±0.06 

13.06 2.16 ±0.11 

13.23 2.62 ±0.10 

13.33 2.95 ±0.12 



Hell 4686 75.50 3.15 ±0.14 
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Fig. 5. — Spatially resolved spectra for the sections defined as "disk" (90° - 270°, light grey) 
and "gas stream" (0° - 90°, dark grey) for a set of six annular regions (labels in units of 
Rli)- Dashed lines show the best-fit HI model for the disk spectra. Vertical dotted lines 
mark the major transition lines of the spectra. 
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Fig. 6. — Top: Radial intensity distributions of selected optical lines (filled symbols) and of 
the corresponding adjacent continuum (open symbols). Dashed lines indicate the la limit 
in each case. Second row from the top: Net line emission radial dist ributions. Dotted lines 



Marsh et al. 



(119901 ) . Third from top: EW 



show the I oc r -L5 radial dependence, derived by 
as a function of radius. Bottom: FWHM of the lines as a function of radius. Dashed lines 
show the law v oc r -1 / 2 expected for gas rotating in Keplerian orbits around a white dwarf 
of mass Mi = 0.6 M & . The predicted spin velocity of material channelled by the white dwarf 
magnetic field is too high to appear in the velocity range plotted in the panels (see text). 
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Fig. 7. — The velocity resolved spectra for the four most prominent emission lines in the DQ 
Her spectra comparing the total emission (black) and the asymmetric component emission 
(grey). Vertical dotted lines mark the line center in each case 
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Fig. 8. — Spectrum of the uneclipsed component. Top: total contribution. Bottom: frac- 
tional contribution. The notation is similar to that of Fig. |5j 
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Fig. 9. — Top: The radial run of the brightness temperature of the best-fit HI model to the 
disk spectra (d = 400 pc and (3 = 2°, solid line with filled triangles) and for d = 525 pc and 
(3 = 2° (open circles). The results for two alternative, extreme cases of low (d = 525 pc and 
/3 = 0°) and high (d = 400 pc and (3 = 3.5°) solid angles (see text) are shown as the upper and 
lower dot-dashed curves, respectively. Dotted lines show steady-state disk models for mass 
accretion rates of M — 10 -8 ' , 10 -8 ' 5 and 10~ 90 M yr _1 . The dashed line shows the critical 
temperature above which the gas should remain in a steady, high-M regime, according to 
the disk instability model. Bottom: the radial run of the surface density. The notation is 
the same as for the upper panel. 
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radial distribution falls much more rapidly with radius than the Balmer lines, with an 
a = 3.15, and does not fit the linear relation of the Balmer lines. This is not surprising, as 
significantly different physical conditions are required to produce the H I and He II lines. 



Smakl (119691 ) and iHuand (119721 ) showed that the double-peaked profile of emission 
lines in accretion disks can be accounted for by making two basics assumption: (i) the 
geometry is flat, i.e., the disk emis sion is restricted to the orbital plane, and (ii) the disk is 



in Keplerian motion. 



Smakl (Il98ll ) showed that the line profile is sensitive to the line radial 



intensity distribution and that the a value could be recovered 



profile. From observational evidences (e.g., 



Fabian et al 



jy fitting a model to the line 



19791 : 



Young &: Schneider 



1980 ) he 



found that the emission lines profiles could be represented by power-laws I oar a , where a 
is expected to be in the range 1 — 2. 



Smakl (Il98ll ) gi ves a list of results where a could be estimated for Helium and Bal mer 



series lin es (Z Cha, 



HT Cas, 



Fabian et al. 



Young et al. 



1979 



19811 : DQ Her, 



SS Cyg, 



Cowley et al. 



Hutchings et al. 



1979 



1980 



Stover et al. 



1980 



Dmitrienko fc Cherepashchuk 



19801 ) . In all cases the value obtained for a is in the range 1.7 — 2.2 suggesting that there 



are no systematic differences between objects of different classes of cataclysmic variables. 



Similar values of a are obtained by fitting the line profile o 



Z Cha and the intermediate polar DQ Her. 



Hutchings et al 



the dwarf novae U Gem or 



(119791 ) found that in DQ Her 



the H/3 line has a = 2.2 while Hell 4686 has a = 1.9. These values are inconsistent with 
ou r results. The discrepa ncy may be explained by the low spectral resolution of the data 



of 



Hutchings et al 



(119791 ). which obtained a from line profiles of low resolution (dispersion 
of 48 A mm -1 ) while we extracted a from spatially resolved spectra obtained from data of 
significantly higher spectral resolution (2 A pixel -1 ). 



Marshl (119871 ) analysed the Balmer lines in the dwarf nova Z Cha. In his study the 



differences between the line profiles in the average spectrum are evident. The line wings 
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become wider from Ha to H5 indicating that the radial line distributions become more 
concentrated for increasing transition energy. In quiescence, Z Cha has a relatively small 
inferred mass transfer rate of M — 10~ 12 M yr^ 1 . Simulations s hows that the LTE model 



is not valid in this case, producing a poor fit to the Balmer lines (jMarsh 



19871 ). 



Marsh! (119881 ) fo und a = 1.8 for the B almer lines in the dwarf nova IP Peg in 



quiescence. Similarly, 



Marsh &: Hornd (119901 ) showed that in IP Peg the radial distribution 



of the Balmer lines have similar slopes when the object is in quiescence but that the radial 
distributions are different for distinct lines when the object is in outburst. The H<5 line 
distribution is steeper than H7, that is steeper than H/3, i.e., the a index of the radial 
profile increases with the transition energy. A similar result was found by 



Marsh et al. 



(119901 ) for the dwarf nova U Gem, with H7 and H/3 showing, 



In contrast, for the white dwarf V2051 Oph, 



r espec tively, a = 1.66 e 1.43. 



Saito fc Baptistal (120061 ) found a = 1.55 for H<5 



and a = 1.78 for H7; a decreases with the transition energy in this dwarf nova. During 
their observations V2051 Oph was in an unusual low brightness state (low mass transfer 
rate). Thus, there are observational evidences that the slope of the radial line distribution 
is different for different emission lines. 

In order to understand the physics behind the line formation in accretion disks it is 
necessary to combine new and improved theory with spatially-resolved observational studies. 
Spatially resolved studies from Doppler tomography and eclipse mapping techniques can 
provide high resolution information about the line distribution in the accretion disk. More 
sophisticated spectral models for accretion disk need to include features such as the vertical 
temperature structure and external sources of heat. 
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5.2. Are we seeing the accretion curtains? 

Here we discuss possible interpretations for the asymmetries seen in the low-velocity 
eclipse maps of the H/3 and He II A4686 lines. These asymmet ries cannot be explained in 



terms of gas stream overflow (e.g. 



Baptista &: Bortoletto 



20041 . see their Fig. 5) as they are 



physically unrelated to the gas stream 



spiral shock waves (e.g.. 



Baptista et al. 



traje ctory, nor can they be interpreted in terms of 



2000 



see their Fig. 2) as the morphology, orientation 
and location are inconsistent with those expected for two-armed tidally- induced spiral 
structures. Moreover, it is also not possible to expla in the observed struc tures in terms of a 



Baptista et al 



1998 



, see their Figs. 3 



front-back asymmetry caused by a flared disk (e.g., 
and 8) because this would affect both line and continuum maps, whereas in DQ Her the 
asymmetries are only seen at specific velocity-resolved eclipse maps of particular emission 
lines. Additionally, the disk opening angle is not large enough to result in measurable 
obscuration or limb-darkening effects (see Section I3"72"j) . The asymmetries observed in the 
H/3 and Hell A4686 line maps seem to be of a different nature, related to something not 
previously seen in eclipse mapping experiments. 

DQ Her is the first intermediate polar being eclipse mapped. It is tempting to interpret 
the structures seen in the low-velocity He II and H/3 maps in terms of the projection onto 
the orbital plane of accretion columns/curtai ns formed ove r the magnetic poles of the 



white dwarf by the channeled material (e.g., 



Ferrario 



1996 



In this scenario, due to 



2 The eclipse mapping procedure has the approximate effect of projecting any emission 
from above the orbital plane back along the inclined line of sight to the point where it pierces 
the disk plane. Thus, emission arising at height Z above the disk is projected a distance 
Z tan i towards the back of the disk and is distorted by t he intrinsic azimuthal smearing 



effect of the eclipse mapping method (IBaptista et al. 



d). 



1995). With our choice of azimuthal 



blur width A<ft (Section 13.21) . any point source in the eclipse maps will appear smeared in 
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different opacities for incoming irradiating photons, and energies needed to ionize atoms 
and to excite electrons to the upper energy levels of each transition, the regions of the 
accretion curtains close to the magnetic poles are prone to reprocess the strong irradiation 
by hard x-rays photons whereas regions farther away receive less energetic photons. Hence, 
high-excitation lines such as He II A4686 trace the emission in the accretion curtains regions 
closer to the magnetic poles, whereas the low-excitation Balmer lines trace emission from 
regions farther away from the x-ray emitting magnetic poles. Lines of different excitation 
frame different parts of the accretion flow and, individually, yield only a partial view of 
the emitting accretion curtains. According to this, a better view of the accretion curtains 
may be obtained by combining the complementary high-excitation Hell map with the 
low-excitation H/3 map. An eclipse map obtained by combining the net line emission 
He II A4686 and H/3 low- velocity eclipse maps is shown in Fig. [TTJ A remarkable structure 
resembling a twisted dipole emitting pattern can be seen near disk center. Although 
this figure encourages one to further pursue this line of reasoning, it poses an immediate 
problem. 

Each eclipse map yields a snapshot of the disk brightness distribution averaged over 
the period covered by the 4 eclipses. Therefore, because it comprises a timescale much 
longer than the 71 s spin period of the white dwarf, changes in brightness distribution 
caused by the varying aspect of the accretion curtains with the white dwarf rotation should 
be completely smeared out in the eclipse map — just as in a long exposure picture of a 
moving target. How could one then obtain a coherent image of fast rotating accretion 
curtains in an eclipse map? The answer can be built in terms of the selection in velocity for 
the accreting material performed by the application of spectral mapping techniques. 

The velocity field for gas in the accretion curtains is a sum of its inflowing component 



azimuth by ~ 20°. 
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Fig. 10. — Distribution of the index a in function of the transition energy for the Hell 4686 
and Balmer series lines. The solid line marks a linear fit for the Balmer lines distribution 
(see Section [5TT]) . 




Fig. 11. — Combined net emission Hell plus H/3 low- velocity eclipse map. The notation is 
similar to that of Fig. [2j 
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(along the magnetic field lines) and a rotational component (due to the white dwarf spin 
period). As seen in Sections 14.21 and 14.31 the velocity field is dominated by the rotational 
component, which is significantly larger than the inflowing component. At most rotation 
phases (i.e., curtain orientations), different parts of the curtains show distinct projected 
velocities along the line of sight and their emission will split onto several different velocity 
bins across the emission line profile. Conversely, a specific velocity-resolved eclipse map 
will contain a mixture of different parts of the accretion curtains at different rotation 
phases. This slicing and mixing effect prevents one from obtaining a coherent image of the 
accretion curtains at all velocity-resolved eclipse maps, with the fortunate exception of the 
low-velocity bin. 

When the curtains are aligned with the binary axis -- and only at this orientation 
- the rotational component will be suppressed (because it becomes orthogonal to the 
line of sight), and the curtains will show only the small inflowing velocity component (see 
Section |4~2]) . contributing to the same particular low- velocity bin. We note that two effects 
contribute to further reduce the observed projected velocity: (i) the inflow velocity must 
decrease after the accretion shock until it reaches zero at the white dwarf surface, and 
(ii) for accretion curtains twisted off the orbital plane, regions close to the magnetic poles 
present the lowest projected velocity along the line of sight. Therefore, looking at a low 
velocity bin is equivalent of framing the accretion curtains always at the same two spin 
phases (where the curtains are oriented at an azimuth roughly aligned with the binary main 
axis), allowing one to obtain a coherent image of them. 

A f urthe r question is why the accretion curtains lag the binary main axis by ~ 30°. 



Hellierl (119991 ) applied Doppler tomography techniques to spin-pulse spectroscopy of 
intermediate polars to find that the optical pulsation arises from azimuthally-extended 
accretion curtains, located typically a few white dwarf radii from disk center. His tomogram 
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of the slow rotator PQ Gem shows a twist in the accretion curtains, with material feeding 
onto magnetic field lines ~ 40° ahead of the location at which it hits the white dwarf. This 
is the opposite case of what is seen in DQ Her and the difference may be in the spin period. 
For comparable magnetic field strengths and mass accretion rates a slow rotator (the case of 
PQ Gem, P sp i n = 834 s) picks material at a magnetospheric radius with higher (Keplerian) 
velocity than the accretion curtains. Thus, the accretion curtains tend to step ahead of the 
magnetic pole. On the other hand, in a fast rotator (the case of DQ Her, P sp i n = 71 s) the 
field lines at the magnetosphere are rotating faster than the disk gas being channelled onto 
the field lines (see Section |4T3"|) . In this case the accretion curtains tends to lag the magnetic 
pole. 

6. Summary 

We applied eclipse mapping techniques to spectroscopy covering 4 eclipses of DQ Her 
in order to isolate the emission from different parts of its accretion flow as well as from 
outside the orbital plane. The main results of this study can be summarized as follows: 

1. Velocity-resolved eclipse maps of the H/3 and Hell A4686 lines indicate the 
presence of an azimuthally- and vertically-extended bright spot at disk rim (at 
Rbs = 0.57 ± 0.03 Rli), which is an important source of reprocessing of x-rays from 
the magnetic poles. This is in line with the concl usion that the bright sp ot is the 



main site of the 71 s optical pulsations in DQ Her ( iSaito fc Baptista 



20091 ). 



2. In the inner regions the disk spectrum is flat with no Balmer or Helium lines, 
suggesting that the emission arises from bremsstrahlung radiation by an optically thin 
gas. In the intermediate and outer disk regions the spectrum shows double-peaked 
emission lines typical of a rotating disk gas. The slope of the continuum becomes 
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progressively redder with increasing radius, indicating the existence of a radial 
temperature gradient. 

3. We fit LTE HI emission models to the flux and slope of the continuum in order to 
infer the temperature and the surface density of the emitting gas. The surface density 
increases by six orders of magnitude from the inner to the outer disk regions. The 
temperatures are in the range T ~ 13500 — 5000 K (for R = [0.15 — 0.65] Rli) and 
can be reasonably well described by a steady-state disk with a mass accretion rate 
of M — (2.7 ± 1.0) x 10~ 9 MqHt^ 1 . The strong emission lines with shallow Balmer 
decrement observed in the outer disk regions are in contrast to the low temperatures 
inferred from the slope of the continuum, suggesting an outer disk structure with a 
hot and optically thin chromosphere (responsible for the emission lines) on top of a 
cool, opaque and dense disk photosphere (responsible for the continuum emission). 

4. The narrow and redshifted Call A3934 absorption line in the total light spectra 
and the inverse P-Cygni profiles of the Balmer and Hell A4686 emission lines in 
spectra of the asymmetric component of the eclipse maps indicate radial inflow 
of gas in the inner disk regions (R < 0.3 -Rli), and are best explained in terms 
of magnetically-controlled accretion inwards of the white dwarf magnetosphere in 
DQ Her. We infer projected radial inflow velocities of (200 — 500) km s _1 , significantly 
lower than both the rotational and the free-fall velocities for the corresponding range 
of radii. 

5. The combined net emission Hell plus H/3 low-velocity eclipse map shows a twisted 
dipole emitting pattern near disk center. This is interpreted as being the projection of 
accretion curtains onto the orbital plane at two specific spin phases, as a consequence 
of the selection in velocity provided by the spectral eclipse mapping. This is in line 
with and strengthen the above inference. 
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6. The spectrum of the uneclipsed light is dominated by Balmer and He I lines in emission 
with narrow absorption cores. The line emission is probably from the extended nova 
shell, whereas the narrow absorption likely arises in a collimated and optically thick 
wind from the accretion disk. 

7. A comparison of the line radial distribution for the Balmer lines reveals a linear 
correlation between the slope of the distribution and the transition energy. 
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Reconstruction of FWHM distributions of emission lines 



Here we address the reliability of spectral mapping techniques to reconstruct FWHM 
distributions of emission lines from spatially resolved disk spectra. Our aim is to check 
whether the sub-Keplerian velocities found for the Balmer and He II 4686 lines in DQ Her 
(Section I4.3[) are intrinsic to the variable or an artifact of the mapping method. 

First, we created a synthetic eclipse map of a steady-state opaque accretion disk 
according to the s et of parameters of DQ Her: white dwarf mass of M\ = 0.6 M^ 



(IZhang et al. 



19951 ). white dwarf radius of R WD = 0.0121 R Q ( Marten et all 09951 ) . Roche 
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Fig. 12. — Top: Greyscale plot of the brightness distribution of the steady-state disk used 
in the simulation. Middle: the emission line profile derived from the brightness distribution 
of the top panel assuming a Keplerian velocity field according to the parameters of DQ Her. 
Bottom: The FWHM of the emission line as a function of radius derived from the spectral 
mapping of 23 light curves along the emission line profile of the middle panel. A dotted 
line shows the v oc r~ 1//2 law expected for Keplerian orbits around a white dwarf of mass 
Mi = O.6M . 
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lobe size of Rli = 0.766 R Q ( |Pettersonlll980l ) and an accretion rate of M = 2.7x 10 9 Mp, yr l 



(Sect ion I4.5p . We adopted the geometry of DQ Her {i = 86.5° and q = 0.66, 



Home et al 



19931 ) to simulate the eclipse of this brightness distribution and to generate velocity-resolved 



synthetic light curves for the same set of bins used in the spectral mapping of DQ Her 
(Section 13. II) . In order to construct the light curve of a given velocity bin, we assumed a 
Keplerian velocity field for the synthetic eclipse map and computed, at each orbital phase, 
the flux contribution only for those pixels the projected (Doppler) velocity of which falls 
inside the velocity bin under consideration. The intensity distribution of the synthetic map 
is shown in the upper panel and the corresponding double-peaked emission line profile is 
depicted in the middle panel of Fig. 



Gaussian noise was added to the synthetic light curves to simulate the signal-to-noise 
ratio of the DQ Her data and eclipse mapping techniques were applied to the resulting 
light curves to generate a set of velocity-resolved synthetic eclipse maps (see Section |3T2|) . 
We combined the velocity-resolved maps to derive spatially-resolved spectra for the same 
set of radial annuli used in the analysis of the DQ Her data, and we extracted the FWHM 
distribution of the emission line as a function of radius following the same procedure used 
in Section 14.31 The lower panel of Fig. [12] compares the derived FWHM distribution with 
the v oc r~ 1 / 2 law used to construct the original eclipse map. Uncertainties were derived via 
Monte Carlo simulations with the synthetic light curves. 

The derived FWHM distribution is in good agreement with the Keplerian expectation 
(at the I—a confidence level) at all radii. These simulations show that the FWHM 
distribution of emission lines can be reliably recovered from spatially-resolved spectra of 
accretion disks obtained with eclipse mapping techniques. This exercise indicates that 
the sub-Keplerian velocity distributions observed in the DQ Her lines are intrinsic to the 
variable and not an artifact of the mapping procedure. 
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